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Abstract 

We explore quantum electrodynamics in (1+1) dimensions at finite temperature 
using the method of Discretized Light-Cone Quantisation. The partition function, 
energy and specific heat are computed in the canonical ensemble using the spectrum 
of invariant masses computed with a standard DLCQ numerical routine. In particular, 
the specific heat exhibits a peak which grows as the continuum limit is numerically 
approached. A critical exponent is tentatively extracted. The surprising result is 
that the density of states contains significant finite size artifacts even for a relatively 
high harmonic resolution. These and the other outstanding problems in the present 
calculation are discussed. 
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1 Introduction 



Developing nonperturbative methods of performing computations at finite temperature 
for hadronic systems is soon to become in demand with the advent of experiments such as 
the RHIC to explore the phase transition to and properties of the putative Quark-Gluon 
Plasma. For zero temperature, there is now much activity in also developing Dirac's front 
form |U or Hamiltonian light-cone field theory for studying nonperturbative field theories 
H and extracting, for example, the relativistic bound state spectra. Reviews can be found 
in ||. Since the purpose of this method is to write down a Hamiltonian, diagonalise it 
and thus obtain a spectrum there would appear to be nothing conceptual in stepping 
from there to the computation of thermodynamic quantities: the spectrum, as the basis 
for developing the partition function, has - at least in principle - all the information one 
could require. 

For this reason we tackle in the present work the finite temperature properties of 
quantum electrodynamics in (1+1) dimensions. Over several years and in a number of 
separate works [Q, [5|, 0], the Hamiltonian and spectrum of the theory have been studied 
with Discretized Light-Cone Quantisation (DLCQ). Comparison with what lattice data 
is available J?]] has been encouraging. It is also known that at least some - though not 
all - of the physics of the fermion condensate in the small fermion mass limit || are 
contained in the spectrum even with the trivial vacuum g !• QED(1+1) is thus the ideal 
place to test the light-cone Hamiltonian method for finite temperature and even to make 
some contribution to what little is known about the 'massive Schwinger model' at finite 
temperature. 

The main subtlety in the DLCQ program is the continuum limit. One works at any 
stage in the numerical work with a finite length interval manifested in momentum space by 
a finite total momentum K, or harmonic resolution, available for distribution amongst the 
'partons' of the theory. Numerically, one computes at various K- values and, when feasible, 
extrapolates to extract the continuum limit. For any finite K the resulting spectrum is 
finite and this of course has an impact on thermodynamic quantities in the sum over all 
eigenstates of the Hamiltonian. We study how this works its way into the various ther- 
modynamic quantities such as partition function, energy and specific heat. The partition 
function for various -fT-values varies enormously - however this is not directly a physical 
quantity. We observe that in the energy for low T the results are roughly if-independent 
while above a certain temperature regime the energy becomes linear in the temperature 
with the slope being independent of K. This leads, in the specific heat, to a if -independent 
result except in the neighborhood of T = 0.5 — 1.0, in units of g/y/n, where a peak with 
increasing height for increasing K appears. The growing peak is suggestive of a property 
of the system as we approach an infinite number of degrees of freedom. We tentatively 
accept this to indicate a phase transition, though it is not associated with the chiral con- 
densate as order parameter (which seems to be nonzero for all finite temperatures anyway 
[pX|]). Rather, it appears to be related to the change in the spectrum from a discrete one 
of bound mesonic states to a continuous spectrum of scattering states. 

First, we briefly review the light-cone formalism as applied to QED(1+1), discuss the 
basic spectrum and give the corresponding density of states as determined in DLCQ. 
We work in the canonical ensemble and present results for the partition function, energy 
and specific heat for various values of harmonic resolution. The letter concludes with a 
discussion of the results and the outstanding problems. 
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2 DLCQ for QED(1+1) 



The literature is already very extensive in references to the basic method of DLCQ g, |3| . 
We can therefore afford to be brief. Our light-cone conventions are those of 0]. Light- 
cone quantization involves initialising independent quantum fields at equal light-cone time, 
say x + = (x° + x 3 )/\/2 = 0. The orthogonal x~ = (x° — x 3 )/y/2 is the longitudinal or 
light-cone space. The appropriate time evolution generator is P~, itself a constant of the 
motion under evolution in x + , together with the momentum or ^"-translation generator 
P + . Thus these operators can be built via the energy- momentum tensor from the bare 
Lagrangian fields once all redundant variables have been eliminated. In order to regularise 
the infrared, a finite interval is employed: x~ £ [-L, +L\. Bosons are assigned periodic and 
fermions antiperiodic boundary conditions. In the following we ignore zero modes []. Then 
the light-cone gauge A + = is permissible in which the remaining A~ gauge potentials 
are constrained and can be eliminated by solving the Gauss law constraints. Thus, apart 
from the resulting linear Coulomb potential, there is no trace left of the photons in two 
dimensional QED. Similarly, only half the fermion degrees of freedom are independent, the 
rest being eliminated by the Dirac equation - a strictly light-cone peculiarity. One is thus 
able to completely specify the Poincare generators in terms of the so-called 'right-mover' 
fermion field which is quantised canonically. The vacuum is just the perturbative Fock 
vacuum in this framework |jl2|. Hence states can be built up by application of creation 
operators in the right-movers after an expansion in plane wave modes. 

The relativistic bound-state problem can be formulated as the eigenvalue equation 



2P + p-\*i) = mf^i) . (1) 

Here, P^ is expressed in terms of the Fock operators into which the fields are expanded, 
and \ is similarly built from Fock operators on the trivial vacuum. Inserting a complete 
set of Fock states results in a finite matrix equation 

(i\ :2P+P- :\j)(j\y)=m 2 (i\y), (2) 

which has been diagonalized numerically @, ||. The reader can find the explicit form of 
the Hamiltonian in [Q, ||. We give here just the conventions for scaling of the coupling and 
invariant mass. With g the gauge coupling and rap the fermion mass in the QED(1+1) 
Lagrangian, we define a coupling 



A = 

and a mass scale 



vrm F 

g 2 



2 I" 1 



(3) 



Q 2 

m\ = ml + — . (4) 

7T 

The coupling A allows us to explore the weak (A — > 0) to strong (A — > 1) coupling regimes 
within a finite range plot. In the following, all energies shall be in units of as shall be 
the units of temperature as well. 

On the computer it is more convenient to express the continuum limit in terms of 
the 'harmonic resolution' K = LP + /tt which is related to the total (integer) momentum 
available to be distributed to the partons of the theory. Evidently, in order to take L — > oo 

1 Their complete incorporation into the analogue of these calculations is still under investigation. 
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Figure 1: Mass Spectrum from DLCQ. a) The Mass Spectrum normalized to uie 
against the coupling A as obtained from DLCQ for antiperiodic boundary conditions and 
a harmonic resolution of K=20. b) The same spectrum but over a smaller mass range with 
larger resolution to emphasise the transition from the discrete to continuum spectrum. 



with fixed total momentum one must then take K — > oo. For some physical quantities it is 
feasible to compute for a range of values of K and extrapolate, however in the present work 
we are restricted to studying thermodynamic quantities for finite but increasing K. Herein 
lies one problem to which we will have cause to refer frequently: computing with K ^ oo 
generates an error in the numerical results, however as yet this error lacks a physical 
interpretation. Nonetheless, the maximum value of K chosen enables reproduction of 
results for the low energy mass spectrum consistent with lattice gauge theory but still 
permits computation on a workstation. The programs are described in some detail in g. 



3 Spectrum and Density of States 

Before diving into the thermodynamics, we prepare ourselves by reviewing the basic prop- 
erties of the spectrum of QED(1+1). A typical spectrum can be found in [|3j]. For com- 
pleteness we reproduce one used to calculate the quantities of this letter; applying slightly 
different boundary conditions as discussed in O. 

There are a number of basic features common to all values of coupling constant. The 
lowest bound state is essentially a fermion-antifermion pair which becomes the 'Schwinger 



boson' of mass g/y/n as fermion mass mp vanishes [ 13 1 . For non-zero fermion mass, this 
mesonic state has a 'size' and so we refer to it as an 'extended' Schwinger boson. Above 
this state lie several bound 'molecule' states of the extended Schwinger boson state. The 
number of these depend on the strength of the coupling and the value of the 6 parameter 
which is essentially a background electric field E3. Finally there is a a continuous spectrum 



consisting of extended Schwinger boson states in relative motion. Embedded within this 
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Figure 2: Density of States in Spectrum from DLCQ. The Density of States in the 
mass spectrum as obtained from DLCQ with increasing values of harmonic resolution: (a) 
K=15, (b) K=20, (c) K=25. The coupling is relatively strong A = 0.85. 

continuum are discete bound states of larger 'molecules' of extended Schwinger bosons, 
for example one of three bosons, namely predominately a six-particle state [ [i"5f| . 

Numerically, for any finite K we always obtain a discrete spectrum which in the limit 
should become a continuum. Put another way, for any K ^ oo there are no scattering 
states, rather states of progressively weaker binding as K increases. In particular, the 
spectrum on its own is rather misleading in illustrating how close one is to a continuous 
spectrum. We give the density of states for a particular value of coupling for various values 
of K in Fig.|. 

For other couplings the results are qualitatively similar. One observes that at some 
point in the spectrum the density of states begins to decrease and this is purely an artifact. 
The point at which the density begins decreasing does not appear to vary strongly with 
K within the range of values explored here. It is difficult to draw any conclusions on this 
aspect in the K = oo case. We do learn here that we should not trust the thermodynamics 
for temperatures commensurate with energies above the point where the density of states 
decreases. This will be reflected in the results to be given in the next section, where we 
take a rather conservative range of temperatures. 

These features aside, we conclude this section by observing the basic structure of the 
density of states: it is basically a rapidly rising distribution breaking down at a certain 
temperature, thus roughly resembling a Gaussian with height and width which increase 
with increasing K. This will be useful at the end in order to understand certain properties 
of the thermodynamics to which we turn now. 

4 Canonical Ensemble 

As mentioned, the DLCQ method involves a calculation of the spectrum for a fixed 'size' 
of system, which is subsequently varied in order to approach the continuum. One should 
note that this size is not related to that measured in any conventional frame accessible 
by proper Lorentz transformations, and so it is not intuitive to picture this as being in 
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a certain frame of reference. Rather, the method tends to a system of infinite length in 
a standard reference frame. Put another way, the results of DLCQ for, say, a spectrum 
acquire relevance to the 'real world' with no boundary walls only when there is some 
converging behaviour as K — > oo. 

Now, as we attempt to describe the thermodynamics of such a system we are lead 
to a canonical ensemble of such finite size systems, there being no external source of 
particles in this problem. Again, these canonical ensembles bear no resemblance to any 
finite size system in a meaningful rest frame. Rather, by taking K — > oo we (hope to) 
approach the same thermodynamics of a set of canonical ensembles of progressively larger 
volumes as determined in the usual sense. Thus, it is wrong to interpret the results for 
any finite K. Rather we must always (where possible) extrapolate at least qualitatively 
to the continuum. 

We now compute the partition function in the canonical ensemble using the basic rela- 
tion Z = J2 i exp(—f3Ei) where Ei is the energy measured in some well-defined rest-frame. 
We consider the system in some given frame attainable by proper Lorentz transformations, 
and now relate the energy Ei to the relativistic-invariant mass squared of the particles in 
the system. Here we make our first genuine approximation: we shall treat all the states 
over which the partition is sum is taken as free discrete bound states, enabling us to take 
M 2 independent of P + and to use the continuum energy-momentum relation: 



Z= ]T £ e-V^F^ . (5) 

states i momenta pi 

The momenta pi are discrete given that the system is taken to be of finite size d ^ L. 
They are the momenta of bound states measured in the given frame for which the size is 
d. Eq.(|5|) is valid for the low energy spectrum which is indeed that of free bound states, 
either of a single extended Schwinger boson, or the molecule states. We now use for Mf the 
values obtained from a DLCQ calculation for a given coupling A and harmonic resolution 
K. Insofar as for any finite K all the DLCQ states are discrete, Eq.([|) would appear to 
be correct. Where the error lies is in concluding from finite K results the thermodynamic 
properties of the finite size system in some Lorentz frame. Unfortunately, at this point we 
cannot determine the error being made following this route. 

From Eq.(|5|) the procedure is now straightforward. A useful trick is to approximate 
the momentum sum by an integral, 

d 



-TMiKWiMt) (6) 
7T t-r 1 



with K\ the first Bessel function in the notation of [16], where we add that integrating 
over P + in the same approximation results of course in the same expression. We take 
Boltzmann's constant to be one so that T = (3 . The energy e and specific heat c are 
calculated via the relations 

dlnZ 2 
e = • T and 

dT 

c = * (T) 
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a) T b) T c) T 

Figure 3: Partition function, Energy and Specific Heat. Results of a DLCQ cal- 
culation with antiperiodic boundary conditions and coupling A = 0.30 for the following 
quantities: a) Partition function, b) Energy, c) Specific heat. 




Now, given the spectrum for a given K and A the above quantities are straightforwardly 
calculated. We first present the results for an intermediate value of coupling and various 
values of K in Figj3|. 

The most significant feature we observe is the rising peak in the specific heat as K 
increases. This region aside, for other temperatures the different K values coincide. Sec- 
ondly, we see a rise in the energy above this peak temperature with increasing K. 

For different couplings, the results qualitatively do not change. Quantitatively, the 
changes are most evident by looking just at the specific heat. In Fig]3] we compare this 
for three different couplings. The only slight change is in the position and height of the 
peak as the coupling increases. More precisely, inspection shows the temperature at which 
the peak occurs appears to approach T = 1 as A — ► 1, which coincides with the strong 



coupling or alternately massless fermion case. We mention that in our units the lowest 
state in this regime, the Schwinger boson, has mass M = 1. 



5 Discussion 

Before attempting to understand these results, we first note that we have given values in 
temperature which do not correspond to systems in the ensemble with states occupying 
deeply into the putative continuum spectrum. Insofar as, with increasing K, we trust the 
low energy derived spectrum the only error we could be making is in treating even the low 
part of the scattering spectrum on the same footing as the discrete one. As mentioned, 
we are unable at this point to quantify how significant this error is. 

Can we nonetheless understand the origin of the peak? Based on our earlier observation 
that the density of states roughly resembles a Gaussian distribution, we can try simple 
analytic calculations using 

p(E)=ex V (-(E-E ) 2 /A) . (8) 

This indeed gives a qualitatively similarly increasing peak in the specific heat with de- 
creasing A. In other words, the height of the peak comes from the steepness of the rise in 
the density of states as K increases. 

The question now is: can these results be extrapolated to the continuum to draw 
physically relevant information? This is impeded by several features. The first is a con- 
ceptual problem: relating if to a physical length scale for a finite size scaling analysis. 
The second impediment is entirely a practical one: going beyond the present values of K 
while maintaining quick CPU time. The question is whether there is some convergence to 
a finite value for the specific heat at the turning point, or whether indeed it diverges to 
infinity. The associated question is whether the profile of the rise in the density of states 
actually converges to some fixed form in the continuum limit. Either way, increased com- 
puting power is necessary which could be complemented by algorithms such as the Lanczos 
method [jT^| . Only via such computations could we say with any confidence whether the 
peak is indicative of a phase transition. 

Let us, nonetheless, assume scaling behaviour, and use the values of the peak in the 
specific heat to attempt a numerical extrapolation over increasing values of K. In this 
way we have estimated the critical exponent a (see, for example, |l8|]) associated with the 
suggested second-order phase transition. Direct fits ranging from K = 10 to 25 resulted 
in increasing exponents from a ~ 0.38 to 0.61. At best, if this is critical behaviour we 
estimate a lower bound for the exponent at a > 0.7. A scaling analysis using K as our 
relevant length scale (in the absence of any other choice) resulted in a determination of 
a and T c for A = 0.05 to a = 0.89 ± 0.04 and T c = 0.54 ± 0.04 consistent with the above 
bounds. Of course, all this is severely limited by the problems discussed above but we put 
these estimates forward as the basis for further discussion. 

At any rate, for low to moderate temperatures the method of calculation would appear 
to be consistent and yields physically sensible results. An interesting calculation which 
could be pursued with the present method is that of the temperature dependence of 
the chiral condensate. There exist analytic results in the literature which allow some 
comparison [10]. 



In summary, we have presented a method for extracting thermodynamic quantities 
from the spectrum of a given field theory as computed using Discretised Light-Cone Quan- 
tisation. The particular example we chose was QED in (1+1) dimensions. Were the DLCQ 
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program for QCD to be achieved, namely the generation of a hadron spectrum, then the 
method presented here could be generalized with little effort. Of course, computing the 
spectrum of QCD from a 'first principles' calculation remains the difficult challenging task, 
but nonetheless the extension of the method proposed in the present work will be essential 
for such computations. 
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